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Finite-well potential in the 3D nonlinear Schrodinger equation: 
Application to Bose-Einstein condensation 
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Abstract. Using variational and numerical solutions we show that stationary negative-energy localized 
(normalizable) bound states can appear in the three-dimensional nonlinear Schrodinger equation with a 
finite square-well potential for a range of nonlinearity parameters. Below a critical attractive nonlinearity, 
the system becomes unstable and experiences collapse. Above a limiting repulsive nonlinearity, the system 
becomes highly repulsive and cannot be bound. The system also allows nonnormalizable states of infinite 
norm at positive energies in the continuum. The normalizable negative-energy bound states could be 
created in BECs and studied in the laboratory with present knowhow. 

PACS. 45.05.+X General theory of classical mechanics of discrete systems - 05.45.-a Nonlinear dynam- 
ics and chaos - 03.75.Hh Static properties of condensates; thermodynamical, statistical, and structural 
properties 



' 1 Introduction 

■ The nonlinear Schrodinger (NLS) equation with cubic or 
. Kerr nonlinearity appears in many areas of physics and 

■ mathematics [T]. Of these, two areas have drawn much 
I attention in recent time. They are pulse propagation in 
, nonlinear medium [SlElSj and Bose-Einstein condensation 

' (BEG) in confining traps [S] . A quantum-mechanical mean- 
, field description of BEG is done using the nonlinear Gross- 

■ Pitaevskii (GP) equation which is essentially the NLS 
' equation with cubic nonlinearity. Though the NLS equa- 
tion in these two areas have similar mathematical struc- 

' ture, the interpretation of the different terms in it is quite 
I distinct. In BEG it is an extension of the Schrodinger equa- 

■ tion to include a nonlinear interaction among bosons. In 
' optics it describes electromagnetic pulse propagation in a 

nonlinear medium. Also, usually there is no external po- 
tential in the NLS equation in optics [1] , whereas in BEG a 
trapping potential is to be included in it [51 . In most stud- 
ies in BEG an infinite parabolic harmonic potential has 
been included in the NLS equation which simulates the 
infinite or nearly infinite experimental parabolic magnetic 
trap. 

In this paper we consider a finite square- well trapping 
potential in the NLS equation with cubic nonlinearity. Al- 
though, we consider a square-well potential for obvious 
analytical knowledge about this potential, most of our re- 
sults should be valid for any finite potential and the ex- 
periments are really carried out on finite traps. This po- 
tential is piecewise constant and leads to analytic solution 
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in many one-dimensional (ID) problems, and serves as a 
model for a trap of finite depth and can be realized in the 
laboratory. Because of these features, there have been sev- 
eral studies of the ID NLS equation with a finite square- 
well and other simple potentials. Zakharov and Shabat [3] 
found the solutions to the NLS equation with a constant 
potential on the infinite line. Garr et al. solved the NLS 
equation under periodic and box boundary conditions [6] 
as well as with the ID square- well potential JJ. There have 
been studies of the step function , point- like impurity 
potential [1D| and the parabolic potential [TT] in the NLS 
equation, as well as transmission of matter wave across 
various potentials [T^ . 

We extend the above ID investigations to the spherically- 
symmetric three-dimensional (3D) square- well interaction. 
This is possibly the most-studied problem in linear quan- 
tum mechanics and allows analytic or quasi-analytic treat- 
ment in many cases. Also, it models an experimental situa- 
tion which can be realized with present-day BEG technol- 
ogy with the use of a detuned laser beam of finite intensity 
[Till] . Such an optical device could generate a square- well 
potential in one direction [7] . Three such potentials in or- 
thogonal directions could make an excellent model for a 
finite 3D square-well potential. In BEG three standing- 
wave orthogonal laser beams have already been used to 
make a 3D periodic optical-lattice potential [13] . In view of 
this, the creation of a 3D square-well potential seems pos- 
sible. Once a BEG is materialized in a square-well poten- 
tial, it could be studied in the laboratory and the results 
compared with the prediction of the theoretical models 
based on the GP equation, thus providing stringent tests 
for these models. 
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We show that it is possible to have normahzable sta- 
tionary BEC bound states in locahzed finite 3D square- 
well potentials for a range of nonlinearity parameters. A 
too strongly attractive nonlinearity parameter is found to 
lead to collapse, whereas a very strong nonlinear repulsion 
does not bind the system. In addition to the normalizable 
stationary bound states, the repulsive NLS equation with 
square-well interaction is also found to yield nonnormal- 
izable states where the probability density has a central 
peaking on a constant background extending to infinity. 
Obviously, these nonnormalizable states do not satisfy the 
boundary condition that the wave function iplr) at a ra- 
dial distance r should tend to as r ^ oo. The forma- 
tion and the study of the normalizable states could be of 
utmost interest in several areas, e.g., optics 1 , nonlin- 
ear physics and BEC [5j , whereas the nonnormalizable 
states will draw the attention of researchers in mathemat- 
ical and nonlinear physics. We use both variational as well 
as numerical solutions of the NLS equation in our study. 

In this connection we mention that in an exponentially 
decaying finite potential well one could have the interest- 
ing possibility of quantum tunneling and the appearance 
of quasi-bound states, which has been studied in detail in 
Ref. [14]. In the present study with square- well potential 
this possibility is not of concern. 

In Sec. II we present the theoretical model which we 
use in our investigation. In Sec. Ill we explain how to ob- 
tain numerically the usual normalizable solutions of the 
NLS equation with the finite and infinite square-well po- 
tentials. We also explain the origin of the nonnormaliz- 
able solutions and how to obtain them numerically. Then 
we develop a time-dependent variational method for the 
study of this problem. The nonlinear problem is reduced 
to an effective potential well. The possibility of the ap- 
pearance of stable bound states in this effective potential 
for a wide range of the parameters is discussed. In Sec. IV 
we consider the complete numerical solution of the NLS 
equation for a finite and infinite square- well potentials. 
We obtain the condition of stability of these bound states 
numerically and find their wave functions. We also ob- 
tain the nonnormalizable solutions of the NLS equation 
numerically. Finally, in Sec. V a brief summary is given. 



where G = At:Ti^ aN / m is the nonlinearity, and v{TZ) is the 
square-well potential. Here m is the mass of each atom, 
N the number of atoms and a is the atomic scattering 
length. We introduce convenient dimensionless variables 
by r = n/l t = TLu /2, ^ = <?>Z3/^ V{r) = 2v{H)/{huj), 
where I = yjh/{muj) and where w is an external reference 
angular frequency. In terms of these new variables the GP 
equation becomes 



d 



Q^2 



2d_ 

r dr 



+ V{r)+g\^{r-t)\^ 



where g = SiraN/l. The square well potential is taken as 

V{r) = -72, r<A; = 0, r > yl, (3) 

with 7^ the depth and A the range. Equation ^ with 
cubic nonlinearity is the usual NLS equation often used 
in problems of optics and nonlinear physics and will be 
referred to as the NLS equation in the following. If we set 
g — Q 'm Eq. ([2]), this equation becomes the usual linear 
Schrodinger equation. In BEC t denotes time and r the 
space variable. In nonlinear optics, t denotes the direction 
of propagation of pulse, r denotes the transverse direc- 
tions, and W refers to components of electromagnetic field. 
In nonlinear optics the 3D NLS equation is spatiotemporal 
in nature where for anomalos dispersion the time variable 
can be combined with the two space variables in transverse 
directions to define the 3D vector r with r = |r|. There 
have been many numerical studies of the 3D NLS equa- 
tion in optical pulse propagation [4[|16|. In BEC a scaled 
nonlinearity n is often defined by n = g/{%TT) — Na/l. 
The normalization condition in Eq. ^ is 



dV|tf'(r;i)p = 1. 



(4) 



3 Analytic Consideration 
3.1 Normalizable Solution 



2 The Nonlinear Schrodinger Equation 

We begin with the radially-symmetric time-dependent qu- 
antum-mechanical GP equation used to describe a BEC 
at K [5] . As we shall not be concerned with a particular 
experimental system, we write the GP equation in dimen- 
sionless variables. The radial part of the 3D spherically- 
symmetric GP equation for the Bose-Einstcin condensate 
wave function r) at position TZ and time r can be 
written as [TS] 



n'dR. 



v{n) 



G\Hn-T)\' 



<^'(7^;r) =0, 



(1) 



The localized normalizable solutions to nonlinear equation 
([2]) with potential Q are allowed only at negative ener- 
gies. Wc solve numerically Eq. ^ starting from a time 
iteration of the linear problem obtained by setting 5 = 
in this equation. Hence we present a brief summary of the 
linear problem in the following [17 . The stationary solu- 
tion of the nonlinear equation ([2]), which we look for, has 
the form W(r]t) = tplr) exp{^iiit) with fi the chemical 
potential, so that 



2d_ 

r dr 



V{r)+gmr)\' 



ipir) ^ fj.ip{r). (5) 



For the piecewise constant square-well potential, the 
solution of the linear problem with g — is expressed in 
terms of the variable a — ^ (7^ — for r < A\ and 
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P = ^y\JI\, for r > A. In terms of these variables the 
solution of Eq. ([5|) regular at the origin is expressed as 



-ip{r) = A 



sin(Q;r) 



--5-cxp(- 
pr 



r < A, 



-pr), 



> A. 



(6) 
(7) 



These solutions are discrete and normalizable satisfying 
Eq. (ID) corresponding to a negative ^. The unknown pa- 
rameter /i is obtained by matching the reduced wave func- 
tion riplr) and its derivative at r = yl: 



aAcot{aA) = -pA. 



(8) 



Once ^ is known the wave function ?/'('') can also be deter- 
mined by implementing the normalization condition ^ of 
the wave function on Eqs. ^ and 

The boundary condition ([5]) is simplified for an infinite 
square-well potential, where sin(ayl) = 0, so that a — 
jTr/A, where the integer j = 1 corresponds to the ground 
state, and j > 1 to the excited soliton states. The solution 
in this case is given by Eq. ([6]) with ip{r) = for r > A. 
After obtaining the solution of the linear equation with 
g = in this fashion, the nonlinear equation ^ is solved 
by time iteration. 



3.2 Nonnormalizable Solution 

The nonnormalizable solutions to Eq. ([5]) with the square- 
well potential ([3]) are only allowed for a repulsive system 
for positive ^ values. The stationary bound-state solution 
of the linear problem discussed above is normalizable. The 
solution of the nonlinear equation generated from that so- 
lution by time iteration is also normalizable. However, for 
positive g (repulsive system) Eq. ([5]) also has nonnormal- 
izable solutions with no counterpart in the linear problem. 
They are obtained from time iteration of a special nonnor- 
malizable solution of Eq. ([5]) for V{r) — 0. The normal- 
ization integral ([4]) is now infinite even for a system with 
finite number of particles TV with a finite scaled nonlin- 
earity n. In other words a system with a finite n (finite 
scattering length a and finite number of atoms N) can 
possess a wave function with nonzero probability density 
everywhere in space. 

We note that 4'{r) = c, a constant, is a solution of Eq. 
([5]) with V{r) = 0, provided that ^ — gc^. It is realized 
that for a repulsive system fj, is positive and the present 
state is a state in the continuum. The required solution 
of Eq. ^ for a nonzero V(r) is obtained from this solu- 
tion by time iteration while in each step of time iteration 
the strength 7^ of the square-well potential is increased 
slowly until the desired value of 7^ is attained and the 
final wave function is obtained. The final wave function 
tends towards the nonzero constant c for r 00. 

For a large condensate the kinetic energy term in Eq. 
(O is negligible and one has the following Thomas-Fermi 
(TF) approximation to the wave function 



As V{r) is piecewise constant, in the TF approximation 
7/;'Pp(r) will also be piecewise constant. Usually, to imple- 
ment the TF approximation we need the parameter 
But in this case we determine /Lt by requiring that in the 
TF approximation V'TF (^) ~^ asymptotically. The actual 
numerical solution also tends to this asymptotic limit. 

If the solution of Eq. ^ is generated from the initial 
constant solution ip{r) = c, the TF approximation to the 
wave function becomes 



^TfW = - V(a^+W3, r<A, (10) 
ilj-[p{r) ^c = ^/Jjjg, r > A, (11) 



ith 



C 



1 7^ + c^g 



(12) 



3.3 Variational Result 

To understand how the normalizable bound states are 
formed in the square-well potential we employ a varia- 
tional method with the following Gaussian wave function 
for the solution of Eq. ([2]) [TBI 



ip{r, t) = A{t) exp 



ia{t) 



(13) 



where A{t) = [■K^''^R^/'^{t)]-'^ , R{t), P{t), and a{t) are the 
normalization, width, chirp, and phase of ip, respectively. 
The Lagrangian density for generating Eq. ^ is P^FIITB] 



V{rm''-\g\^\ 



dip 



dr 



(14) 



where the overhead dot denotes time derivative. The trial 
wave function p3p is substituted in the Lagrangian den- 
sity and the effective Lagrangian calculated by inte- 
gration: LgQ- = / C{'ip)(Pr. The Euler-Lagrange equations 
for this effective lagrangian are [l9] 



d ^-^eff ^ ^-^eff 
dt dq{t) dq{t) ' 



(15) 



where q{t) stands for the generalized displacements R{t), 
P{t), A{t) or a{t). 

The following expression for the effective Lagrangian 
can be calculated after some straightforward algebra [18] 



^eff - 



71^/2^2 (t)i?3(t) 



2a{t) - 



3P^{t)R^{t) - TAjiR) 



(16) 



where 
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with the error function Erf(a;) defined by 

2 /"^ 2 

Erf(x) = ^ / e'*^dt. 



(18) 



The Euler-Lagrange equations (|15p for a{t), A{t),(3{t), 
and R{t) are given, respectively, by 

^3/2^2^3 ^ constant = 1, (19) 

2gA' 



4q, 5 2 



(20) 



R = 2RP, (21) 



16 7^^^ 



30F i?5 



r/r^-.VH^ = _4^, (22) 



V2i?2 



where the time dependence of different observables is sup- 
pressed. Ehminating a between ((20|) and ((22|) one obtains 



16 -f^A^ 



(23) 



From pTjl and (|23p we get the following second-order dif- 
ferential equation for the evolution of the width R 



fR 



dU{R) 
dR 

4 8n 1 
i?3 ^ 



i?4 30F i?4 



^'^\-A^/-\ (24) 



where n = g/{8TT) and ?7(-R) is the effective potential of 
motion given by 



R^ 3V2^ i?3 3 



(25) 



Small oscillation around a stable configuration is possi- 
ble when there is a minimum in this effective potential 
denoting a stationary normalizable state. 

In Figs. 1 (a) and (b) we plot U{R) vs R of Eq. 1^ in 
dimensionless units for different values of n, A and 7. We 
exhibit U{R) vs R for different n for yl = 7 = 2 in Fig. 1 
(a). The same for different 7 and A = 2 and n = 10 are 
exhibited in Fig. 1 (b). We find from Fig. 1 (a) that for a 
fixed 7 and A, U{R) has a minimum for n^^^ > n > n^j-it 
corresponding to a stable bound state, where n-^[jj^ is pos- 
itive (repulsive system) and n^^.^^ is negative (attractive 
system). In Fig. 1 (a) there is a minimum for n = 10,0 
and —0.4 and no minimum for n = 40 and —1. The reason 
for the nonexistence of bound state for n > n^^jj^ and for 
n < 'T.fjj.j^ are distinct. The limiting condition n > n^^^ 
corresponds to a highly repulsive system which cannot be 
bound in the square-well potential with the given A and 
7. The condition n < n^j-j^ corresponds to a highly at- 
tractive system which collapses with the given A and 7. 
In Fig. 1 (a) for n = 40 the system is unbound and for 
n = — 1 it is unstable to collapse. In Fig. 1 (a) we see that 
as the attractive nonlinearity \n\ is increased, one of the 




Fig. 1. The effective potential U{R) of Eq. ^ vs. R in di- 
mensionless units for yl = 2, (a) 7 = 2 and n = 40, 10, 0, —0.4, 
and -1 (upper to lower curves) and for (b) n = 10 and 7 = 1,2, 
and 3 (upper to lower curves). 



walls of the well for small r is gradually lowered and for a 
sufficiently large attraction this wall is completely absent 
and the condensate collapses into the infinitely deep well 
at the center for rt = — 1 in Fig. 1 (a). 

In Fig. 1 (b) we illustrate the effect of increasing the 
strength of the square-well potential for a fixed nonlin- 
earity n = 10. This is done by varying the depth of the 
square- well potential 7^ for a fixed range A. For 7 = 2 and 
3, U{R) has a minimum and the attraction of the square- 
well potential is sufficient to bind the repulsive system 
with n = 10. However, for 7 = 1 there is no minimum in 
U{R) and the square-well potential is too weak to bind 
this system. 



4 Numerical Result 

We solve NLS equation ^ numerically for the square- 
well potential using the split-step time-iteration method 
employing the Crank-Nicholson discretization scheme de- 
scribed recently '20,21]. We calculated the solutions with 
real-time propagation. However, we checked that imaginary- 
time propagation also leads to consistent result. To obtain 
the normalizable solution, the real time iteration is started 
with the known solution of the linear problem with scaled 
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2 4 6 8 10 
(b) A (|im) 

Fig. 2. Critical nonlinearity n^j.^^ (a) vs. 7^^ for A = 2 fim 
and (b) vs. A for 7^ = 4Eo = Ahuj: full line - numerical; dotted 
line - variational. 



nonlinearity n — Q. Then during time iteration the non- 
linearity n is switched on slowly until its desired value is 
attained. The change in the parameter should be such that 
it does not greatly alter the eigenvalue of the Hamiltonian 
(after time propagation). We also calculated the nonnor- 
malizable bound states in the continuum for a positive n. 
To obtain this solution, the time iteration is started with 
a constant wave function for a finite positive n with 7 = 0. 
Then in the course of time iteration the strength 7 of the 
square-well potential is switched on slowly until its de- 
sired value is attained. If stabilization upon time iteration 
could be achieved for the chosen parameters one already 
obtains the required nonlinear bound state in the square- 
well potential. In previous studies we compared critically 
the present numerical scheme for the time-dependent NLS 
equation with other approaches [20l[2T| including the time- 
independent schemes [22j and assured that the present 
approach leads to results with high precision not only for 
the NLS equation with one space variable but also for NLS 
equations in two and three space variables [23j . The agree- 
ment between the results obtained with real and imagi- 
nary time propagation also assures the correctness of our 
results. 

Although we calculate our results in dimensionless units, 
typical parameters for an experimental realization can be 
easily obtained therefrom for a particular atom. In the 
following we present results for the Rb atom. For Rb let 



>, 180 




5 10 15 20 25 

/(Eo) 

Fig. 3. Limiting nonlinearity ny^-^ vs. 7^ for A = 2 fim. 

the length / be 1 /im; for this to be possible the reference 
frequency is cj w 27r x 38 Hz. Consequently, the unit of 
energy is Eq = Tiuj~ 1.57 x 10~^^ eV. 

4.1 Normalizable States 

Stable normalizable bound states are indeed found in all 
cases for various ranges of parameters. Some plausible 
properties of these bound state are found in agreement 
with the above variational study. For a given nonlinear- 
ity, these bound states are only formed for a sufficiently 
strong square-well potential determined by A and 7. For 
weaker potentials, from the wisdom obtained in varia- 
tional calculation, the effective potential U{R) does not 
have a minimum and there cannot be any bound state. 
For a given square-well potential, bound states are found 
for njipj >n> n^^-^^. 

It is difficult to obtain the limiting nonlinearity n = 
njj-^ numerically as this corresponds to a state with zero 
energy which extends to a very large r. However, the crit- 
ical value be obtained numerically in a 
controlled fashion as the wave function is highly localized 
in this limit. In Figs. 2 we plot the critical nonlinearity for 
collapse for different parameters of the square-well 
potential. In Fig. 2 (a) we plot JT-crit for yl = 2 jim 
whereas in Fig. 2 (b) we plot n^^j.^^ vs. A for 7^ = AEq. 
In these figures we show results of variational and full 
numerical calculations. The variational calculation always 
leads to a larger |n(,j.j^|. In case of the infinite harmonic 
potential also, the variational estimate of is larger 
than the result of full numerical calculation For the 
infinite square-well potential with 7^ = 00 and yl = 2 
/im, ^^CTit ~ —0.62; for the infinite parabolic potential 
"crit ^ —0.575 [5ll24]. Because of the different shapes of 
these two infinite potentials JT-^j-j^ is different in these two 
cases. 

Although it is difficult to obtain ny^^ from a numeri- 
cal solution of the NLS equation, it is possible to obtain it 
from the variational calculation. The limiting nonlinear- 
ity 'T'ljjjj corresponds to the largest value of n for which 
U{R) has a minimum. In Fig. 3 we plot ny^^ vs. 7^ for 
A = 2 /im. We see that ny^j^ increases linearly with the 
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10 



Fig. 4. (a) Ground-state wave function ^li{r) of the NLS equa- 
tion (HI with the square- well potential for 7^ = AEq, A — 2 
/xm and scaled nonlinearity n — 0.75, and 10 (upper to lower 
curves). In all cases the numerical and variational results are 
shown, in addition, for n = 10 the TF approximation is also 
shown, (b) Same for the first excited soliton (numerical calcu- 
lation only) with = 4i?o, A = 4 and n — —0.65 and 6. 
(c) Same for the second excited soliton (numerical calculation 
only) with 7^ = AEq, A = 4 /im and n = —0.4 and 6. 



strength of the square well 7^. The variational calcula- 
tion underbinds the system. The repulsive nonlinearity 
destroys binding and a smaller repulsive nonlinearity can 
destroy the weaker binding of the variational model. Con- 
sequently, the variational limiting nonlinearity is smaller 
than the actual limiting nonlinearity, which we verified in 
our numerical calculation. The numerical calculation re- 
lies on the existence of a localized wave function and it 
is difficult to calculate the limit when this wave function 
extends to infinity and a localized wave function ceases to 
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Fig. 5. (a) Ground-state wave function ■(/'(r) of the NLS equa- 
tion ([5]) for the infinite square- well potential of range A — 2 /xm 
and scaled nonlinearity n = —0.62, and 10 (upper to lower 
curves); (b) same for the first excited soliton for n = —0.3,0 
and 10. 



exist. The TF wave function ^ is always fully localized 
within r < A and is inadequate for calculating the limit- 
ing nonlinearity for /i — > 0. The TF wave function leads to 
a much smaller value n^iYa ~ l^i based on imposing 
[i — Q within the TF regime. 

In Fig. 4 (a) we plot the wave function for the bound 
state of the NLS equation for a square- well potential with 
7^ — 4£'o, A = 2 /im and for nonlinearities n = —0.75,0 
and 10. The system with n = —0.75 is attractive, n = 
noninteracting, and n = 10 repulsive. In these cases results 
for both the full numerical calculation and variational ap- 
proximation are shown. In addition, for n = 10 the TF 
approximation pop with c = /j, = is also shown. In Fig. 
4 (b) we show the wave functions for the first excited soli- 
ton with a single node for 7^ = AEq and A = 4 /itm and for 
n = —0.65 and 6. In Fig. 4 (c) we show the wave functions 
for the second excited soliton with two nodes for 7^ — 4Eo 
and A — 4 /im and for n — —0.4 and 6. In Figs. 4 we find 
that the bound state for the attractive system with nega- 
tive n values is more centrally peaked than the bound state 
for the repulsive system with positive n values. This is true 
for both ground and excited solitonic states in the finite 
square-well potential as well as for states in the infinite 
square- well potential studied below. The central peaking 
of the wave function for the attractive system correspond- 
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Fig. 6. Wave functions of the nonnormalizable states tp{r) of 
the NLS equation ((5]) for (a) 7^ — 4Eo, A — 2 fim and different 
n values, and for (b) n = 0.3, A = 2 and different 7. In both 
cases the wave function tends to 0.1 as r — > 00. In all cases the 
results of the full numerical calculation and TF approximation 
are shown. 

ing to a large central probability density is a consequence 
of the nonlinear attraction. 

Next we consider the infinite square-well potential of 
range A. In this case the system remains confined in the 
region < r < A. The wave function is zero outside this 
region: r > A. We illustrate the wave functions in this case 
for different values of nonlinearity n for the ground state 
and the first excited soliton for yl = 2 /im. For the ground 
state there is no node of the wave function for < r < A. 
For the jth excited soliton these are j nodes of the wave 
function in this region. The wave functions for the ground 
state and the first excited soliton for different n are shown 
in Figs. 5 (a) and (b), respectively. 



4.2 Nonnormalizable States 

Now we discuss the nonnormalizable solutions of Eq. ([2|) 
obtained for repulsive condensates. No such states are 
found for attractive condensates with negative n. To ob- 
tain these solutions by time iteration of Eq. ([2]) we start 
with the initial constant solution ■0(r) = 0.1 over all space 
for V{r) — and a fixed nonlinearity n. In the course 
of time iteration the square-well potential is slowly in- 
troduced without altering the nonlinearity until the full 
square-well potential is obtained. The resulting wave func- 
tion of this calculation is plotted in Figs. 6 (a) and (b) 



together with the TF approximation. In Fig. 6 (a) we 
show the wave function for 7^ = 4£'oj A = 2 ^m and 
n — 0.1,0.3, and 10. In Fig. 6 (b) we show the same for 
A = 2 fim, n = 0.3 and 7^ = Eq, AEq and I&Eq. From Fig. 
6 (a) we find that for a fixed 7^ and A the TF approxima- 
tion becomes a better approximation as the nonlinearity 
increases; whereas from Fig. 6 (b) we find that for a fixed 
A and n the TF approximation becomes a better one as 
the strength of the potential 7^ increases. 

The normalizable states discussed in the last subsec- 
tion are true bound states for negative fi. The nonnor- 
malizable states occur at positive n and can be referred 
to as states in the continuum. Similar nonnormalizable 
states have been obtained in the study of the ID NLS 
equation with the square- well potential [7]. However, there 
is the possibility of physically observing the nonnormal- 
izable states with infinite norm as a transition from the 
normalizable states with very large nonlinearity [7] . An in- 
crease of the scattering length a via a Feshbach resonance 
[25] in a BEC may increase the nonlinearity n(= Na/l) 
indefinitely and thus transform a normalizable state to 
a nonnormalizable one. This conclusion results from the 
following observation. As the normalization integral ([4]) 
and the nonlinearity g — Sim of Eq. ([2]) are mixed in the 
present nonlinear model, an increase in g can either be 
implemented my increasing the nonlinear term g directly 
in Eq. ([2]) or by augmenting the normalization integral (|4]) 
in the same proportion. Hence the nonnormalizable states 
with very large (infinite) normalization could be a transi- 
tion from normalizable states with large nonlinearity. 



5 Discussion and Conclusion 

The study of the ID NLS equation with the simple square- 
well potential is of interest because of its simplicity and 
intrinsic nonlinear nature. Its interest in BEC and optics 
has motivated recent investigation of this problem [7] . In 
the 3D world, ID systems can only be achieved in some ap- 
proximation and there are nontrivial differences between 
the solutions of the NLS in ID and 3D l]. 

Hence we performed in this paper an investigation 
of the 3D NLS equation with the square-well potential 
using numerical and variational solutions. We find that 
the system allows normalizable bound-state solutions for 
""lim > > '^crit' where the limiting nonlinearity n-^^jj^ 
corresponds to a repulsive (positive) limit beyond which 
the system cannot be bound and the critical nonlinearity 
n^jiil to an attractive (negative) nonlinearity below which 
the system collapses. We calculated n^^^^ for different pa- 
rameters of the square- well potential. Many results of this 
paper can be verified experimentally in BEC, where one 
can make a square-well trap by joint magnetic and opti- 
cal control [Zlll]. In addition to the discrete normalizable 
states we find that the 3D NLS equation with the square- 
well potential also sustains nonnormalizable states in the 
continuum of interest in nonlinear physics and mathemat- 
ics. These states cannot be obtained via a transition from 
the normalizable states. 
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